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A brief summary of the literature pertaining to the theory of 


cueues is presented to acquaint the reader with the mathematical 
techniques available. 

A quantitative analysis of the landing process around an air 
craft carrier is made with a view towards a more efficient utiliza- 
tion of the landing platform. The theory of queues supplies the 
mathematical techniques, thus making it possible to come to grips 
analytically with the problem of reducing the landing time. 

The problem has been idealized severely but is applicable under 
certain conditions to the actual landing process. The results may 
be used for investigating the distribution of the means of (a) the 
aircraft waiting time, (b) the time wasted by the aircraft carrier 
during the landing operations, and (c) the number of aircraft waiting 


provided we know (er assume we know): 





(1) the minimum safe landing interval, 

(ii) the distribution of the landing times (landing rate), 

(111) the distribution of the tim intervals between aircraft 

entering the landing pattern (departure rate), 
If we have not assumed the above, this information must be based on 
data obtained from observation of actual Landing operations. 

A method is also indicated whereby the prohlem may be analyzed 
by a quasi-emperical technique. This is a well known substitute for 
an analytical treatment if the mathematics become intractable. 

ii 








With the advent of the jet age, terminal congestion amd resulting 


aircraft delays are of increasing importance in all branches of avia- 
tion. The scientific approach provides a means for predicting how 
these delays will be affected by varying conditions, and so mkes possi- 
ble an economic valuation of efforts directed towards reducing these 
delays. 

Operational research methods have a threefold application in this 
connection: (a) in ascertaining the wide field of practical problems 
susceptible to treatment; (b) in obtaining the basic data needed to 
enable the delay problems to be specified; (c) in selecting approp- 
riate solutions, 

Oueueing theory supplies certain mathematical techniques whereby 
the operations analyst can come to grips analytícally with the problem 
of aircraft delays. Even where mathematics become intractable, there 


are good prospects of empirical solutions by means of electronic 





calculators. 

The purpose of this paper is to apply queueing theory to a prob- 
lem of sans importance in naval aviation, carrier aircraft landing de- 
lays, It follows that the results may allow a more efficient utili- 
gation of the landing platform. Two things are of particular interest 
in the general queueing problem, the size of the queue (stack) and the 
waiting time of the customer (aircraft delay). 

Chapter 1 gives an analysis of the general queueing problem and 
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integral to the problem, throw considerable light on the report. 


ter 11 presents a detailed analysis of one aspect of the carrier air- 


craft problem, a summary of the results, the conclusions, nd further 
considerations of statistical implications, 

This study was undertaken at the United States Naval lostgraduate 
School during the latter hlf of the academic year 1953-54. “he writer 
is deeply indebted to Frofessor Charles C. Torrance whose constant 
criticism and encouragement were invaluable. To him go my foremost 
thanks, The writer is also indebted to Trofessor A. Boyd ilewborn for 


reading the paper and for making valuable suggestions and remarks. 
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CHAPTER I 


QUEUEING THEORY 


1.1 Introduction to the Cueueing Problem, 

The study of queuss has been of interest to mathematicians and 
engineers for the past forty yesrs. ^ueueing occurs vhen demands for 
service are greater than availability of service, and it becomes neces- 
sary to postpone the demands by a system of queueing or mrshalling. 
“hile queueing is a natural process, marshalling is an attempt to do 
something about it. “ueueing, however, is essentially a temporary 
phenomenon, as otherwise the queue would grow indefinitely in size, 

Thus one of the first problems in any queueing situation is to determine 
whether or not the demand will outstrip the service mechanism. Vhen 

this does not happen, so that the queue is continually returning to 

zero size, then the situation is said to be in equilibrium, or, in statis- 
tical terminology, the process is stationary. ‘hen the demand eauals 

the service capabilities the situation becomes very delicately balanced; 

a slight reduction in demand yields a sensible queueing system in equi- 
librium, while a slight increase results in an ever lengthening queue. 

The theory of queues has a special appeal for the scientist interes- 
ted in stochastic processes, for under certain conditions it provides 
an example which is both stationary and non-Markovian. Thus one may 
expect in studying the cueue to gain insight into other stochastic 
phenomena, and to accutire a valuable facility in the handling of the 
relevant technicues, In addition, the theory has an astonishing range 

1 








of applications. A current limitation of the theory is that the num 


ber of papers that deal with the theoretical side of the subject is 
small. And unfortunately those that exist are rather heavy reading. 
1.2 Cueusing system. 

The simple cueueing system is concerned not so mech with the 
matching of supply and demnd, but with a series of customers demand- 
ing service at a single counter and waiting in turn to be served. (us- 
tomer is here used in a technical sense; in the present paper it should 
be ecuated with aircraft. To complete the specification (Kendall 2 ) one 
mist give a careful aceount of: 

(a) the input-process, 

(b) the nueue-discipline, 

(c) the service-mechanism. 
1.2.1 Input-process, 

The input is random whenever the chance or probabillty that a 
customer will arrive at any instant of time is the same for all instants 
past, present, and future, and also is unaffected by any past or future 
arrivals. The random input is often called a "poisson input" because 
of the use of the Foisson distribution to describe the arrivals. A general 
independent input (Lindley 3 ) is one where the intervals between the 
arrivals of successive customers are independent with identical pro 
bability distributions, but this distribution is arbitrary. The input 
is regular when the customers arrive at a constant rate. The input is 
correlated whenever the intervals between arrivals are not independent, 
1.2.2 “ueue-discipline. 

The cueue-discipline is the rule or moral code determining the 

2 








manner in which the customers form up into a queue, and the manner in 





which they behave while waiting. In the simplest case they line up be- 
fore the single counter and await their respective turns. The theory 
of multiple queues or many servers seems, except under simplifying as- 
sumptions which do not always correspond to reality, to be a problem of 
considerable difficulty and will not be considered here. 

1.2.3 Service-mechanism, 

The service-mechanism may be described as the mode of service. In 
this connection the service time is of basic importance and is the time 
that elapses while a particular customer is being served. 

1.3 Traffic density. 

A parameter of special importance is the traffic intensity or den- 
sity € , defined as the ratio of the average arrival rate to the average 
service rate. This may also be expressed as the ratio of the mean ser- 
vice time to the mean arrival interval. 

l.l Markovian properties. 

Before the mathematical theory is presented, an interesting ouali- 
tative feature of the problem should be commented on. 

Suppose that the state of a stochastic system at time t is described 
by a random function X(t); the stochastic process is then said to be of 
Merkov type (Doob 9) if a knowledge of the present value of X(t) makes 
all information about its past history irrelevant to a prediction of its 
future behaviour. Such a process may become brian if part of 
the information contained in X(t) is suppressed, for a knowledge of the 


previous behaviour may enable some of the suppressed information about 
ጋ 














the present state to be recovered. The oueueing process with a Fois- 
sonian input and a single server is Markovian if the present state of 
the process is described by the pair of random variables q and vo where 
q is the instantaneous cueue size and YS is the expended service time 

of the customer at the head of the queue; in general it ceases to be 
Markovian if the state of the process is measured by the cueue size 
alone (Kendall 2). The only exception to this statement occurs when 
the service time has a negative exponential distribution. 

1.5 Poisson input. 

Kendall ( 2 ) gives a general review of some points in congestion 
theory. A presentation is made of the Pollaczek "equilibrium" theory 
for the single counter queue fed by an input of the Poisson type and 
associated with a general service time distribution. It is pointed out 
that although the stochastic processes describing the fluctuations in 
queue size is not in general Markovian, it is possible to work instead 
with an enumerable Markov chain if attention is directed to the epochs 
at which individual customers depart. The ergodic properties of the 
chain are investigated with the aid of Feller's ( 8 ) theory of re- 
current events, 

The simplest hypothesis about the input is one which states that 
the customers arrive "at random", the number of arrivals in time T 
being a Poisson variable of expectation T/a, where a is the mean 
arrival interval, The time interval t between two consecutive arrivals 
will then have the negative exponential distribution 


me 
)عل‎ 2 6 dt 0 > 6 > © , ው 
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This is seen if one recalls that dF(t) is simply the probability of 
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tatistically independent. This is the hypothesis adopted by Kendall. 
1.5.1 Mean waiting time with general service time distribution. 

Let the queue be of the type with a single server fed by a Iois- 
sonian input and let the traffic intensity € be less than unity, so 
that the system is not saturated. Let 4 be the size of the cueue which 
the departing customer leaves behind hin, not including himself, but 
including the next person to be served (a may be zero). The next per- 
son to be served has a service time v [z(v) = b] , Say, and suppose 
that during this time r new customers arrive. “hen conditionally r is 
a Foisson variable of mean value v/a Es 7 ከ88 ህር service time dis- 
tribution. Let ^! be the size of the oueue which the next person leaves 


behind him. 





Assume that statistical equilibrium exists. Then the random var- 
fables q and a! have the same marginal distribution, and in particular 


their mean and mean square values must be eoual and finite: 
E(S):E(S 1 «o, ፡ (ፍና) ፡ 5(ና")< =. 125 
The variables a 5 q' are related by the fundamental formula 
qy = m ax 1 ዘ 1.3 


or in terms of the § notation, 
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tant to note the following consequences of the definition of the 


y. مس‎ UK | 
function ፅ ( qy): 


Ana (A-0) ፡ ay. 1.5 


On forming the expectations of both sides of Eq. 1.4, it will be 
found that 
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This is the chance that $ is non-zero; that is, it is the probability 


that an incoming customer will not have to wait. On squaring both sides 





of Ea. 1.4 and making use of Eq. 1.5, one finds that 


ሙ fon G A \ 
Liar ተ 27-ع -) نه پر‎ 


q * z ዓ - 2 Y (1-0 ) % سوه‎ i (25-1) : እ የሻ 
31|-ዩሥነጆረዬነ 5 A rla Ela + ሯ ሀ”) + ነ፤- ١+ د - نمع‎ 60 
Take expectations and note that r is independent of q and §. It 
follows that 1 
nen 
ጀ(ፍ) ፡ ድ (ጾ) + E نا‎ ] 


2 ١١-٤٤٤ | E. 


በዜ p vanr v æ b‏ د 
e 20 (u-b)‏ 
Now suppose that a departing customer leaves q customers behind‏ 
him, and let his own waiting time and service time be respectively w and‏ 


y. Then a is the number of arrivals in a total time w + v. Thus, 


AL 


Using Eos. 1.8 and 1,9 one finds that 
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This is essentially Kendall's result. It is convenient to express 
the result in the form of the two ratios, for in applying the results to 
aircraft delays, one usually normalizes time in units of the mean arrival 
interval or mean service time (minimum safe landing interval). 

"hen the mean arrival interval and the mean service time are equal 
to each other (€ » 1), the expected waiting time is infinite, thus show 
ing that the intuitive solution € = 1 is not the best result. 

From Ecs. 1.10 or 1.11 one can see that the minimm expected wait- 
ing time will be obtained if the service time distribution is of the 
form de ስብ 


Fiv)=0 (¥< b) ; Fiv)=1 (yb) | 1.12 


that is, the service time is constant. On the other hand, if the service 


time is of the form 
-W 
d F(v) = € dv ١ O< V <% 145 
b 


the expected waiting time is twice the minimum value (var > = 1), 


Tt is tempting to think that the expressions for E(q) and E(w) may 





be valid for more general input processes; but this seems to be a fallacy. 
The difficulty is that r and q are no longer statistically independent. 
7 














Kendall shows that the Markov chain associated with the queueing 
process is irreducible and the states are aperiodic. A further classifi- 
cation of the states depends on the value of the traffic intensity, € . 
"hen € 1s less than unity, all states are ergodic; when € is greater 
than unity, all states are transient; in the critical case when € is 
oqual to unity, the states are all recurrent and null. (Feller 9 ) 
1.5.2 Aircraft delays. 

Bowen and Fearcey ( 6 ) give an elementary analysis of the effect 
of control procedures on the flow of air traffic, while Fearcey ( 7 ) 
presents a more rigorous treatment of the problem and derives the dis- 
tribution of delays under different traffic conditions. The hypothesis 
of random arrivals is used, so that one is led to the negative exponen- 
tial distribution for the arrival intervals. See q. l.l. 

The assumption is made that the minimm safe landing interval 
(service time) is constant, Thus we have Kendall's solution with a con- 
stant service time in a different notation. 

Fearcey's treatment is based on a rather complicated procedure of 
successive integrations but once the assumption is made that the system 
is in a state of statistical ecuilibrium the analysis mst lead to the 
same results. Cne feature of interest is the derivation of the proba- 
bility of successive delays. By successive integrations Fearcey finds 
that in general the probability of n-l aircraft being successively 


delayed is (stack of n-l aircraft) 
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Bell ( 5 ) gives a somewhat more general treatment by using a moment 
generating procedure. 
1.6 General independent input, 

The center of interest of Lindley's ( 3 ) paper is the waiting times 





of the customers, 1 relation between the wiiting times of successive 
customers is found, and this makes it possible to calculate the wait- 
ing tim distribution of any customer, 

This relation is also used to obtain necessary and sufficient con- 
ditions for the existence of a stationary (e-uilibriun) state. The paper 
concludes with an investigation into the stationary state in some spec- 
ial cases; particularly in the case when the customers arrive at regular 
intervals. 

1.6.1 Waiting time distribution. 

The „ueue here is again of the type with a single server but with 
general service time and input distributions; the intervals between 
successive arrivals my have any suitable distribution. 


Let t. be the tire interval between the arrival of the a and 


h customer. 


r + ¡e customer, and let 5. be the service tire of the m 
The following assumptions are made: 
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v^ 


l. The b. are independent random variables with identical 
probability distributions and the mean, E(t..), is finite: 


ES. [e em ete, - 1 


2. The 3 are independent random variables with identical 
probability distributions and the mean, =(8..)› is finite, 
Thus the two sets of random variables EN and CV] (P 1,2,7 
are statistically independent, 
Let w, be the waiting time of the r'! customer. It is evident 
that the following relations exist (Fig. 1.1): 


P. | n ቀ! — Y annıv | 
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Fig. 1.1 ‘Waiting Time Relations 
دم م۷ + با‎ ۲ ۷ ተ 9. (tp < m * 8L) 1.1, 
t 
i.a: T (tr 2 m + 5,5 see t.) 1.15 


Let «دپا - پوه په‎ The u are independent random variables 
(by virtue of the proceding assumptions) and have identical probability 
distributions; E( | u, | ) is finite and u, is independent of w,. Then 


* 


et (m * U, > 0), 1.16‏ و با 
(m, * U, <0). 1447‏ 0 = 
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Now if we assume that the first customer does not have to wait 


(w, = 0) we have the following relations: 


i ۲ WOES Uu E ۱ 
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and the unexpected fact develops that the waiting times depend on only 
the difference between the service time and interval time, If we define 
F(x) as the probability that w, $ x, and G, (x) as the probability that 
انا‎ <x, we can slightly modify Lindley's procedure and find the distri- 
bution of waiting times of the r + ገኙ customer by using Ea. 1.18. The 
u are independent statistically. Thus the distribution function of 


r 
F, + l is the convolution of هېه لا‎ tee tu (Doob 9) and 
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Using the definition of G.(x), we have 
FA. (x): p (u rutet un SX : ا‎ UTK, rsh) 1.20 


Eq. 1.19 establishes a recursive relation such that the individ- 
ual waiting time distributions my be found. To find the final behavior 
of the queue, we mst find the limit of CO as the number of cus- 


tomers continually increase. Let E, be the event: 


ifs 
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U + u,+ u U, t + USE vee 1.21 


Now the sequence of events ES decreases to a limit E as r approaches 
infinity, that is: 
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by a theorem of probability measure. In the limit, then, 


Fix) = | F(x-u) dG (u) 1.23 
ما‎ 5 
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This integral equation is difficult to solve, but it is possible to do 
so once G(u) is specified and certain conditions are satisfied (Paley and 
Weiner ll‘). The solution of either integral equation gives the dis- 
tribution function of the waiting time, that is, the probability of 
waiting a time not more than x, The distribution function for any par- 
ticular customer can be calculated by performing the necessary inte- 
grations; but this is also difficult to do. 

Lindley uses Feller's theory of recurrent events and the strong 
law of large numbers to investigate the behavior of F(x). 

The event of "not having to wait", called event A, is a recurrent 
‚event, since the queueing process may be considered to make a fresh start 
whenever a customer does not have to wait, and its previous history is 
quite irrelevant to subsequent developments, 

To investigate the behavior of F(x) one mst examine the relation 


ሮ (ሠነ ። ፒ (5) = E («&) 1.25 


as follows: 
E(u) > 0: The strong law shows that the event A is transient and 
12 











that with probability 1 there will be som point in the queueing process 


where every subsequent customer will have to wait. 

E(u) <0: The strong law shows that event A is certain with a finite 
mean recurrence time. It should be clear that it is not periodic. 

E(u) = 0: The event A is certain, but the mean recurrence time is 
infinite. Hence it is certain that there will always be another occasion 
on which a customer does not have to wait but the probability of any 
particular customer not having to wait tends to zero. 

Thus a necessary and sufficient condition that the waiting tim dis- 
tribution function tends to a non-degenerate limit, as the number of cus- 
tomers increases, is that E(u) < O or u = O certainly. | 


Lindley salves the integral equation when the arrivals are random 


that is, 
. - ۸ سف‎ 
(ې)6 4ل‎ «© dry, Ey)’, 1.26 
and also for 
dG6(4): Aye "ta, Ely): 2 
Y m T 55 ) Y X 1027 


These distributions are the Type 111 curves of Pearson commonly called 


the garma distributions (Mood 10), and are of the form 








N+! በ 
: ^ 07 4 = - ገግጎተ 1 1.23 
with the given parameters n and À . 
And as expected, the probability of not having to wait is 
TERN) ع داد‎ 
۱ E(t) 1:23 
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The equation to be solved for regular arrivals, where the arrival 


interval is the unit on the time scale, is 





Fl): Jug F Orgel de, (yer): 1.30 
with the service time distribution as given by Eq. 1.28, 
Nti x N = ኣኣ 
Faa fry et 
) m! 0 4 4 ك‎ : 1.32 
The solution is difficult to develop, but many may be expressed in the 
form 
Nel ۳ xX 
ም (።=ንነ 5 |- GRON 
۰ ብ ) 1.32 
ዬኒ = = and 
ገጎ 
كد‎ e E SE 2:0 رت او عم‎ 
a ; +1 4 J 1 ) 
5 (አተይ ) 
ዶድ:። £. and 
NU" ሺ a- Z 
nat | e 


The usual form of service time distribution is one whose frecuency 
function is zero at 3 = O, increases to à maximum and then decreases, 
with a long tail, to zero as s approaches infinity. In some cases it 
may decrease steadily from some non-zero value at s = O, To a good 
approximation these functions are represented by the garma distribu- 
tions (n = 0, 1). 

The results of the calculations show that regular arrivals are 
roughly twice as efficient as random arrivals, i.e., the probability of 
not having to wait for regular arrivals is roughly twice as great as 
the random arrivals, Lindley's calculations for a few values of € are 
listed in Table 1 to show the nature of the results (units are arbitrary). 

1), 














Mean waiting time 





ال دج rt‏ 8 پا په € 
4.166 1:823 0.167 0.314 0.830 
1.786 ,0.681 0.286 0.51 0.715 
1.042 0.349 0.375 0,642 0.635 
0,500 0.128 0.500 0.797 0,500 
1 = ۲ 
ጋር‏ 0.834 0.167 0.386 0,830 
1:999 0.285 0.286 0.608 0.715 
0.781 0.132 0.375 0.743 0,635 
0.375 0.039 0.500 ,0.884 0.500 
Table 1‏ 


1.7 Correlated input. 

In view of the diffieulties encountered with the random and inde- 
pendent inputs, it should be readily apparent that the inclusion of ae 
ditional probabilities in the input process will make the problem even 
more difficult. As might be expected then, this problem has not been 
solved analytically. 

The use of the Poisson input in air traffic leads to a complete 
analytical solution of the stack-delay problem, but unfortunately once 
aircraft are in the air, the input is not sven approximtely randon, 
but has the form of a regular pattern with large errors. It is virtually 
impossible to find the actual distribution with any accuracy, as the 
amount of data required for an unbiassed sample is prohibitive. The key 

| K 








to the difficulty lies in the fact that aircraft in the air proceed on a 


predetermined schedule. This detailed information in the schedule should 
be regarded as a datum of the problem and used in its solution. 

The first significant contribution of the electronic calculator to 
the queueing problem has been its use in the study of the flow of 
scheduled air traffic. The Research Laboratory of Electronics at the 
Massachusetts Institute of Technology ( 4 ) are the pioneers in this 
endeavor. They undertook a program to determine, as far as practicable, 
the quantitative relation between the degree of control of aircraft 
enroute to an airport, and the resulting congestion and delay. 

Actually there are two possible problems involved in the control 
aspects of air traffic. One is the maintenance of a pre-assigned schedule, 
and the other is the reduction in the congestion around the terminal or 
landing strip. It should be noted that this is an unusual situation in 
that the aims of the customer and the server are the same: (a) to ensure 
safety in the air; (b) to allow the most efficient utilization of the 
landing platform; (c) to reduce the total time by avoiding long delays 
at the terminal, This is not always true in the general aueueing problem, 

In the MIT report, IBM punched-card machines were employed to make 
a theoretical analysis of the resulting congestion at a single landing 
strip when aircraft are scheduled to arrive there in some proper ee- 
quence, but fail to meet such schedules according to certain simple 
deviation statistics (rectangular, triangular, and parabolic). These 
deviations statistics are assumed to have a finite spread, that is, 
deviations in excess of a specific amount have zero probability. This 
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is within the bounds of reason from the pilot's viewpoint. Numerical 


results for the distribution of the resulting stack and total delays 
are presented. The former are compared with those which arise from 
the more random Poisson arrival distribution. The amounts of conges- 
tion and stack delay are found to be considerably smaller under the 
new conditions, particularly when the traffic is heavy. A significant 
range of the parameters has been covered by the numerical analysis, 
supported by strictly analytical methods wherever possible. Thus there 
is established a quantitative relation between the time-keeping errors 
of the aircraft enroute to an airport and the resulting terminal con- 
gestion and delay. 
| The density of the traffic flow to an airport is described by the 
usual traffic parameter €, defined as the ratio of the average arrival 
rate at the airport (number of planes per unit time) to the maximum 
allowable acceptance rate. This is also the ratio of the landing in- 
terval (service time) to the arrival interval. 

The authors found that it was convenient to use time in discrete 
units rather than on a continuous basis. The minimum safe landing 
interval, t,, is chosen as the unit of time, and forms the base for 
both the normalization and quantization of all other times appearing 
in the problem. 

A sample size of 1,000 planes was selected for use in the numerical 
analysis. In the programming procedure, each plane in the sample is 
represented by a card, which is identified by a number j, with j greater 
than or equal to l and less than or equal to 1000. Thus j defines the 
scheduled arrival order. The scheduled arrival times t 3 are then as- 

m" 








signed using a random number table. 


Between its scheduled take-off and actual arrival times the aircraft 


ኔክ card ) 


is subject to a delay rj, which is assigned (punched onto the j 
according to the enroute distribution used, Thus plane j actually ar- 
which number is also entered on the jth 


rives at time Pj = t, + r 


t.‏ ال 
card, It is possible now for more than one plane to arrive at the‏ 
sams time and for them to arrive in a se-uence which differs from that‏ 
scheduled. This leads to the formation of a stack of aircraft, and‏ 
consequent delays while waiting in the stack. It is then assumed that‏ 
planes will be handled at the terminal in order of actual arrival, with‏ 
the provision that simultaneous arrivals will be considered as ordered‏ 
according to their originally scheduled sequence.‏ 

After the preliminary details were completed, i.e., the t ۲ and r 1 
tables were prepared, and the initial conditions were chosen, the IHM 
machines were programmed in detail for the specific cases, 1.e., for 
each deviation distribution, a spread S, and traffic parameter € . 

The principal results which came directly from the machines were: 

(a) Freruency distributions for the stack delays. 

(b) Average stack delay for each run. 

(c) Cumulative frequency distributions for the stack delays. 

(d) Fre-uency distributions for the total timekeeping errors 

(stack delay plus enroute delay). 

(e) Average timskeeping error for each run. 

(f) Cumulative fre-uency distributions for the total timekeeping 

error. 
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Without loss of generality, it was assumed that all variations 


were delays (simple change of scale). 

The results form a considerable source of data pertinent to a 
traffic of properly scheduled aircraft. The data are given in detail 
in the form of curves of stack-delay probabilities, prcgressive-delay 
probabilities, and average delays, for the three types of enroute-de- 
viation distributions considered. 

The curves give additional information regarding the effects of 
variation of the parameters. In particular, only relatively small 
differences in the results are caused by rather moderate variations 
in the shape of the enroute-deviation distributions. The relatively 
small. changes in the stack-delay distribution caused by making moderate 
changes in the shape of the enroute-deviation for given S made it 
possible to average together the TEM results for all three basic shapes. 
The average curves consequently give stack-delay probabilities for an 
average shape enroute-deviation distribution, which is sort of cross 
between the rectangular, triangular, amd parabolic shapes. 

The effects of changes in the traffic parameter € are practically 
the same for all three types of distribution. For values of € close 
to unity, the average delay changes fairly rapidly with € , while cor- 
responding changes in € , when it has smaller values, cause little 
change in the average delay. The stack-delay distribution is also de- 
pendent in & similar way upon the traffic parameter. ‘hen € is close to 
unity, the most probable stack delay is almost equal to its average value, 
but as € decreases, the most probable delay rapidly approaches zero 
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with a small average delay, The remaining variable parameter, the spread 


S, also has a decreasing effect with decreasing values of € 
1.8 Tabulation of cases covered. 

The principal cases covered may be tabulated briefly as follows: 

(a) Poisson input with general service time distribution (Kendall 2 ). 
Derives expressions for the expected waiting time and the ex- 
pected queue size, 

(b) General Independent Input i.e., general service time and input 
distributions ineluding regular arrivals (Lindley 3 ). 
Derives the waiting time distribution of any customer and the 
همتا‎ function of the waiting time in the stationary 
case. 

(c) Correlated input i.e., scheduled air traffic (Adler and Fricker 
با‎  ). Use electronic calculator (IBM) to find the frequency 
distributions for the stack delays, cumulative frequency dis- 
tributions for the stack delays, and the average stack delay 
for each run. The information is presented in the form of 
useable graphs, 

1.9 Implications 

An interesting interpretation of the MIT results may be visualized 
in the following manner. Consider the problem of bombing when a large 
number of planes may be directed to strike several targets. The usual 
practice, to avoid saturating the handling capacity of the airfield, is 
to stagger the formations in such a manner that all of the aircraft will 


not return at the same time. 








For example, consider a 200 plane strike made up of 20 groups of 


10 planes per group. Now we have no idea of the enroute deviation dis- 
tribution, but let us use our imagination and say that the average 
enroute-deviation distribution described in Section 1.7 is applicable. 

Assume the expected flight time to target and back to be 10 hours 
for each group. Let the minimum safe landing interval be approximately 
2 minutes per plane so that it takes on the average 20 minutes to land 
a group (t, = 20 minutes). 

It is desired that the 20 groups arrive home in such a manner that 
the planned maximum stack is to be 2 groups or less (20 planes). Larger 
stacks are not feagible due to lack of control facilities, safe altitude 
separation, etc. 

For a 10 hour flight, errors of * 20 minutes would not be deemed 
too large; in the notation of the LOT report, this will correspond to a 
spread of 40 minutes. If we normalize the spread in units of t., S = 2. 

If we apply the results of the MIT paper with € =], that is, with 
the departure interval also 20 mimites and with S = 2, we find the pro- 
bability of stacks greater than or equal to 2 groups is approximately 
0.95. In other words, it is almost certain that more than 2 groups will 
overlap at the field. If we reduce € to 0.8, that is, make the depar- 
ture interval 25 minutes, the probability that stacks greater than or 
equal to 2 groups will form is reduced to 0.06. Thus it may be seen 
that queueing theory, as developed by electronic computers used in con- 
junction with good experimental data, has striking possibilities. 

This example is of course trivial in that it is such a restricted 
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ein. But even so, the results are well worth further study. The 


Poisson solution when € is l has been stressed enough; for example, the 
stack approaches infinity, and the waiting time approaches infinity. 
Thus there would be no solution in that event. hen € is 0.8, the 
Poisson solution gives for the probability of 2 or more groups in the 
stack a value of 0.25 approximately. In comparing the Foisson solu- 
tion with the IBM solution we see the two results differ by a factor 

of h4, a considerable discrepancy. The key to the difficulty is the 

use of schedule information as a datum of the problem. 


$ 
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CHAPTER 11 


CARRISR AIRCRAFT LANDING DELAYS 
2.1 Formulation of the problem., 

The prime desideratum in landing aireraft aboard a carrier is 
simply to land the aircraft as quickly as possible while observing the 
necessary rules of safety. The tims required to land (and launch) the 
aircraft aboard the aircraft carrier is of basic importance due to the 
extreme vulnerability of the carrier during the landing operations. The 
intuitive solution of the problem of minimizing the landing tim is to 
minimize the landing intervals by feeding the aircraft into the landing 
pattern as quickly as possible subject to the restraint of minimm safe 
landing interval. ‘hile this is sometimes a satisfactory solution, yet 
sometimes the delays involved becom excessive. Thus in attempting to 
minimize the landing tins, one must investigate the behavior of the air- 
craft during the landing process. This fact makes necessary a quanti- 
tative analysis of the landing process, with a view towards a more 
efficient utilization of the landing platform. Thus it seems natural 
to apply the theory of ousues in the analysis. Let us consider the 
main aspects of the situation under the three headings introduced by 
Kendall to characterize a queueing process. 

As a preliminary step, consider the pattern of events around an 
aircraft carrier. All of the aircraft waiting to land are placed íin an 
orbit above the carrier at a safe altitude until the signal to begin 
landing is given. This may be visualized as a ready-made stack (queue). 
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"hen the signal to land is given, the aircraft descend to the landing 


pattern to commence their approach. There is a constraint in that only 
a maximum number of aircraft are allowed in the landing pattern. Every 
time an aircraft lands it is replaced in the Landing pattern by one of 
the waiting airoraft. Thus one can visualize the single server queus- 
ing process. 

First, there is the input process. In our context this is essen- 
tially the process by which the aircraft depart the orbiting stack (queue) 
and feed into either the landing pattern directly or into another stack 
(queue) of aircraft waiting to enter the landing process. The landing 
pattern may also be visualized as a queue with its associated pre 
blem of delays dus to wave offs, etc. This is a problem for further 
study and is not considered here. Thus we may visualize two queues 
with a single server. It should be noted that the orbiting stack (high) 
is called a queue in a strictly limited sense. This artifice makes the 
analysis more understandable. The low stack (queue) is the feature of 
interest in that it is caused by some perturbation in the landing process. 
We shall assume that all the aircraft leave the high stack at a con- 
stant rate (continuous flow), i.e., at a succession of equally spaced 
times, the difference between one time and the next being called the 
departure interval. Thus the departure distribution is of the ۳ 


۳ ۴ (۳ 0 (rcy) ۱ Sve (nev) 


where r is the departure rate (constant). Second, there is the queue 
discipline. Ve shall require that the aireraft land in the order in 
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which they arrive at the low stack if one exists, Third, there is the 


service mechanism, which is given by the frequency distribution of the 
landing times (serving time). In a limited sense we shall consider the 
landing pattern as a part of the service mechanism. Landing time is 
here taken as the time from the cut signal given by the landing signal 
officer to the time the aircraft is safely forward of the barriers. ۵٥ 
will now make the drastic simplifying assumption that the landing time 
is constant, i.e., the landing distribution is 


,)4*5( ا: (و) چ ۰ (5 >۸) © Feel):‏ 


where 1 is the landing rate (constant). This assumption is mde for 
two reasons: (a) observational data are lacking, and (b) it is desired 
to make the analysis tractable. It seems rather natural to assume that 
the landing time distribution is one whose frequency function is cut off 
at the minimum safe landing interval, increases to a maximum near some 
average value of landing time, and then decreases to another arbitrary 
cut off point as the landing time increases. Some consideration will 

be given later to certain statistical implications pertaining to the 
gathering of observational data, <A second drastic assumption to be made 
in connection with the landing tim is that under normal conditions 
(when the normal Ea rate is in effect) the aircraft are landed 
aboard without any delay for a certain length of time, and then the 
normal landing rate is changed abruptly, and continues at & new reduced 
constant rate for a period of time T,. The reason for this reduction 
in rate may be poor handling technique aboard the carrier, inexperienced 
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pilots, etc, During this period of tim the low stack begins to form. 


It is now easily seen why we consider two queues in series with a single 
server. 

An important parameter of the problem is the number N of planes 
initially in the high stack. Under normal conditions (when the low stack 
is empty) the time for landing N aircraft is 

፲ጉ = N/r. 
Thus the problem is to investigate the delays under the two conditions, 
> < 
To 2 T. and Ts چ‎ 1 
Before proceeding with the analysis it should be noted that the 


r. #2 


fundamental difference between the general queueing problem and the pro- 
blem at hand is that we are not interested in the limiting case, that 
is, in an infinitely long queueing process. The stack (queue) assoc- 
lated with the carrier problem has à finite bound N as a datum of the 
problem, and this must be included in the analysis. However, as the 
departure interval approaches the landing time, we might expect from 
the work of Lindley that the situation is somewhat critical for short 
imate well as for long ones. This is in fact the cage. | 

The notation to be used is as follows: 

Q(t): Landing rate aboard aircraft carrier (planes per unit of 

time). 
r (t): departure rate from high stack (planes per unit of time). 
€ : normal traffic density parameter also equal to the ratio 


of the landing tim to the departure interval . 


9 


landing rate reduction parameter. 
T : total tims a is diminishing the landing rate. 
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T, : total time during which non-zero low stack exists. 


Tn : total time for carrier to complete landings under normal 

conditions. 

N : number of planes initially in high stack 

S(t): number of planes in low stack (low stack size). 

T__ : total additional time required for carrier to complete 

landings due to formation of lew stack. 

Y : aircraft delay. 

D : total aggregate delay (say in plane-mimutes) equal to product 

of average delay (7 ) and the total mmber of planes delayed; 
also ecual to the product of T, and s(t). 
Note that the traffic density € is not a function of time. 

The numerical quantities of primary interest will be: (1) the ad- 
ditional time required for the carrier to complete landings due to the 
reduction of the landing rate; (2) the average stack delay over the per- 
lod during which a non-zero low stack exists; and (3) the maximum stack 
delay suffered by any plane. 

The calculations fall into two m jor groups, defined by b 3 Ta and 
T$ Tye THe first case to be considered will be T. ፲.. Figs. 2.1, 
2.2, 2.5, and 2.6 show the corresponding departure rate, landing rate, 
and stack height for prescribed values of o and € in the ranges 

> گر > ۶ ۶ 9 

We can assume without any loss of generality that T 9 commences 
when t = O, Prior to this time, the landing process is normal, say 
(€: Y ን €< 1 ) and there is no low stack. It is not essential to 
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the analysis whether the landing process is in progress at t = 0, or 


just commencing. 

Now at the beginning of time Tg the landing rate is reduced to al 
and, provided a «€, a low stack begins to form at a rate (r-a &).. 

r-a = 2(€ -a). 

It is essential to keep in mind that no stack forms at all unless a <e, 
and this condition is therefore implicit throughout the derivations. At 
time وې‎ the departure rate from the high stack becomes zero, and, in 
consequence the stack no longer increases but begins to decrease at a 
rate af . After a further period of time (T ) the aircraft carrier 
resumes its norml landing rate. As the high stack has cleared, the low 
stack decreases with the normal landing rate ۰ 
2.2 Stack sige S(t) and stack time پل‎ 

At time T, the low stack height is 


5 (ኺ.)፡ ኺ ዳዩ -=) 2 243 


S(T,) is ecual to S(T) ax. ' See Fig. 2.1. Now the stack begins to de- 
crease at a rate al, and at time To, 
ST): T Lte-a)- ۲۵۰ 1 ( ۷ ۸ 
> TAR 6-۵ tor, ). "a 


Now S(T.) is non-negative if and only if 


"le ) 
if this is not the case and 


é > T, 
a 


F< lc 
the stack clears entirely prior to the end of time To’ Ses Fig. 2.2. 
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Thus we have two cases to investir 
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The tims taken to clear the stack left at time 2 15 
كد‎ (€-x lo ) 
[ E 
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Now the time that a non-zero low stack exists is 


- T 2 - o le 
Ts o + Tn te "7. 2.5 





ل كن 


The stack height at time T, is 


SC): TR (e-«). a 
In this case the stack begins to decrease at a rate © / ; thus the time 
for the stack to clear is 
Samy) E TER A (c-a) 

R XX 


په 
2.6 


RSA (&--!) 
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Average low stack S. 


iv 


Case A: ር 2 16 
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The average stack height by definition is 
ا‎ Ts 
oa DN 1 S(t) dt 


Ts 
ue | ET 
. LT 1)» اه (ه‎ dt 


* (fa (e-a Te ) - tajat f, 


and after integration and simplification 
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Case B: 


Ey c 5 
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R 
The average stack S is easily seen from Fig. 2.2 to be 
> 2 1 ۳ پچ‎ = 
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2,3 Average delay T . 


The derivation of the average delay by the usual process for com- 
puting averages is rather unwieldy. It is mch more convenient to com- 
pute T from the ratio of the average stack height S(T) to the average 
landing rate Y(t) or the average departure rate r(t), all averages 
being taken over the period T, during which the stack is never zero. 
That T can very generally be computed correctly in this manner my be 
demonstrated in the following way. 


S(t) 


+ trdt E FU. + ست‎ 
Stack Height vs. Time Plot for General Case 
Fig. 2.3 
Consider a low stack commencing at t = O and ending at t = Tg 28 


shom in Fig. 2.3. Assume the variation of S(t) as arbitrary but subject 
to the following restraints: 


S CHO, 2‏ ۰ (8) و 


S(t)>0 fon o<t < Ts. 
2,11 
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In terms of the notation previously defined, it follows from Eos. 2.10 


and 2.11 that 


+ 
S(t): | Incer-Re)] dt 


© 
o r OT 


This of course puta certain general restrictions upon r(t) and ሐ i5), 
since they must lead to & stack subject to Eqs. 2.10 and 2.11. Granting 
these general restrictions, Eq. 2.12 is valid regardless of the detailed 
form of either r(t) or 1 (t). Now since S(T,) = 0, 


at || f^ (t)dt a 


That is, for such a stack variation, the total number of planes depart- 
ing the high stack equals the total number landed, 

Consider the planes r(t)dt which arrive in the time interval be- 
tween t and t + dt. They will have to wait say, until t = t4 to land, 
because of the time required to clear the stack S(t). It fellows there- 
fore that t4 is a function of t established by 

k 
S(t): | Alt) dt 

t 2.1, 
which simply states that all the planes in the stack at time t mst 
have landed at a rate (| (t) by the time tj. From Eqs. 2.12 and 2.14, 
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7 In (€) - git) | dt : | ፈ6] dt 


E IT 
J nit) dt : | Rit) dt 


2.15 


which defines u as a single valued function of t, Thus as might be 
expected, the delay of a plane is a function of the time it arrives. 

Now by definition, the average delay "T 1s the total aggregate 
delay D divided by the number of planes delayed, where 


b: [P (£-t)i(t)de,- 


That is, the planes arriving in the tims interval dt are delayed a time 
(t, - t) and the delays must be summed over the period T, during which 
the stack is never zero. nd of course the number of planes delayed is 
simply the total number arriving during the period Tg» Thus 


i T, ) 
[ ۲۳ ۱ 4+ 2.16 


۱ or in view of Eq, 2.13, 
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PFT 2.17 
J * R(t) dt 





Now the numerator can be integrated by parts as follows: 
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By Eqs. 2.10 ard 2,15, 
t; = O when t =0, 2.19 


ty = Ta when t -T,- 2.20 


This simply states that the aircraft which arrive when the stack is zero 
are not delayed at all. Using Eq. 2.15 to eliminate 
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from the last term of Eq. 2.18, we have 
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and by a similar process 
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Thus 
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Now dt, is merely functioning as a variable of integration in the last 
term of Eq. 2.23, so that 
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Using the rather obvious definitions 
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and Eas. 2.12 and 2.13, we have the desired alternative form of Eos, 
2.16 and 2.17: 
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Using Eqs. 2.4, 2.26, and observing Fig. 2.1, we find that the av- 
erage landing rate is 


E Te $ 
FX" f 1.45 + | Rat 
5 


6 


ER 


2.29 











Thus from Eos. 2.28, 2.29, and 2.8, the average delay is 
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From Fig. 2.2 and Eqs. 2.26, 2.28, and 2.9, the average delay is 
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2.4 Maximum delay ይኔ ۰ 


In order to find the maximum delay, the delay of planes arriving 
at all possible times within T, must first be considered. 


Stack 
Height 
S(t) 


A. 


€ 4 + یرل‎ To Ts = ዴኤዴኤዴኤ-፦- ጐ-፦ 
Fig. 2.4 Stack Height vs, Time 
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As illustrated in Fig. 2.4, a group of planes arriving in the time 


interval (t, t + dt), witht $T_, finds a stack of / (€ - a)t planes 
waiting to land. Now the landing rate is a X at any time t € To, so 
that the delay for the group of planes arriving during time dt is 


m: t & («- <<) . t (&-&) 
eg j ሚ 2,32 


ዕ ፈጭ ሩ We 


Eq. 2.32 is only valid if the delay Y $ Ty = t, for ifT>T,-t, all 
of the planes will not have landed by Tg and we have to consider the 
change in landing rate at T,. Thus another limiting condition on t in 
Eq. 2.32, in addition to the condítion 

6 ۶ ۲ > 7 
is that 
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The real limiting condition on t in Eq. 2.32 therefore depends upon 


whether 
7 T To 2.31, 


A 
Te se, 2.35 











For the condition 


۲. ۶ * ۸ (case A) | 
2.36 


the planes arriving in the interval dt will not have landed by T,. 
So the number of planes still due to land before them at a. is 


ቲ ፪ርዩ-ጪ)- (CT, -t) Kk. 


2.37 
Because these planes land at a rate |, 
Au (T-t) + RES - UTE EY. 
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Now it should be clear that the maximum delay will occur for some 
group of planes arriving st or just before Tr, where the stack takes 
on its maximum value and the landing rate remains small. The difficul- 
ties of finding the maximum near the end point of a linear ecuation 
of the form of Eo. 2,38 need not be elaborated on. The method used 
here is empirical and the value of t obtained to maximize Ec. 2.38 is 


با 
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For the condition 
t iO INS = To (case 8) 
ፎ / 


2.40 


the delay for any plane in the group is given simply by Eo, 2.32, and 
for T to be a maximum set 


t eT 
Thus 
: E 
Tar 7 TA (= 5 
2.41 
for To < ۴ < | 3 Case B, 
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2.5 Additional delay of aircraft carrier Tie 
It should be clear that the extra time needed by the aircraft car- 
rier to land aircraft vhen a lew stack exists is simply 
Tae " La - Te 2.42 
where 
EN N/r. 
Thus T,, 48 the additional time reruired for the carrier to steam into 
the wind at high speeds and should be made a minimum for the most ef- 
ficient utilization of the landing platform. 
Case A: £ > Ts > | 
A Ta 
From “as. 2.4 and 2.42, 
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From Eqs. 2.7 and 2.42, 
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2.6 Summary of results. 

In the same manner, for የር E the calculations split up into ۰ 
two portions, depending upon the time at which the stack finally clears, 
See Figs. 2.5 and 2.6, The complete results are summarized in Table 2, 
which gives the length of the stack period, the maximum stack size, the 
average delay, the maximum delay, and the aircraft carrier delay. 

From Table 2 the aircraft carrier delay Tq, is seen to be in general 


a function of the variables T5, T., a , and &. It is convenient to 


وم 
use a normalizing factor in preparing graphs for visual representation,‏ 
le was chosen as the normalizing factor for convenience because in prac-‏ 
tice it is usually fixed. A fanily of curves, for a fixed value ofa ,‏ 
showing T,./Tr vs the variable y = T/T, (to specify Tg) for‏ 
several values of € as a parameter, were prepared. See Figs. 2.7, 2.8,‏ 
and 2.9, Other curves could be prepared in a similar manner.‏ 

The general implications of the numerical results may best be ap- 
preciated in terms of examples. Consider a carrier conducting landing 
operations with the traffic density € approximately l, Now let us 
assume that during the landing operations the landing rate has been 


suddenly reduced from the normal rate in such a manner that when a = 0.7, 
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other words during landing operations the landing rate has been reduced 


in several instances so that on the average, say 0.8 of the time, land- 
ings were being made in an inefficient manner (poor pilot technique, 
inexperienced plane handlers, etc.). Now from Figs, 2.7, 2.8, and 2,9 
we find that the additional delay to the carrier Tac to be approximately 
OobTr. That is, the carrier must spend almost half the time normally 
taken to land aircraft steaming at high speed using extra fuel ín a 
most inefficient manner, Taken over a period of time this could amount 
to a sizeable wastage, to say nothing of the danger that might be 
present in enemy waters, 

Similarly, other values of the parameters may be used to compute 
new delays, and in addition the behaviour of the stack may be inves- 
tigated, i.e., maximum stack, average delays, etc, These may be impor- 
tant in determining fuel reserves of the aircraft. 

The hazards in working with the traffic density € approximately 
l have been stressed throughout the paper. However, carrier landing 
operations are a unique application of queueing theory and it should be 
clear that the most efficient utilization of the landing platform will 
only be accomplished with € approximately 1, And this can be accomplish- 
ed only when the pilots and other personnel are exercising the most pre- 
cise skill, judgment, etc. during the landing operations. This is of 
course merely a statement of what common sense would dictate. 

2.7 Conclusions, 
For a simplified situation occuring in the landing of carrier air- 


craft, an investigation has been made into the variation of the average 
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aircraft waiting time, the time wasted by the carrier during landing 
operations, and the distribution of the number of aircraft in the stack 
for different initial conditions, The simplifications and assumptions 
necessary to arrive at a simple solution make the problem somewhat arti- 
ficial, However, it is not completely unreasonable, because in actual 
practice the landing process may actually behave in some instances in 

a manner similar to the process assumed, 

The practical verification and application of the results des- 
cribed in this paper depend upon the specification of the landing time 
distributions, the traffic density © , and the minimum safe landing 
interval. 

In comparing the present results with actual carrier operations, 
many observations of this process must be made during actual landings 
so that the validity of the assumptions made may be verified. In prac- 
tice, then, a great deal of information is required to check the 
numerical model which has been used as a basis for the result pre- 
sented here. If these assumptions are not valid as based on observation- 


al data, the model discussed in the next section may be applicable. 
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2.8 Statistical Implications. 


The difficulties encountered in any analytical püti of a queue- 
ing process with regular arrivals and a respectable distribution of ser- 
vice timea is apparent from Lindley's analysis (Sec 1,6). In particular, 
if the process is not in statistical equilibrium, the individual waiting 
time distributions are difficult to obtain (Eq. 1.19). 

An alternative process is a quasi-emperical investigation using an 
electronic calculator as in the VIT report, or if that method is not 
avallable, the use of random numbers may be substituted. This is a well 
known standard alternative to a purely theoretical analysis, and its 
basic principle is as follows. If we read off successive groups of, say, 
four or five numbers from tables of random numbers, then these can be 
regarded as numbers, with four or five decimal places, drawn at random 
from a rectangular frecuency distribution with limits O and 1. Now all 
continuous distribution functions are themselves rectangularly distribu- 
ted between O and 1, so that the abscissa of any given distribution 
function corresponding to a value of the function taken from a table of 
random numbers can be regarded as a cuantity drawn at random, with the 
appropriate frequency, from the frequency distribution itself. Thus for 
a given distribution of landing time, for example a garma distribution 
with given parameters, n and A say, we can use random numbers to run 
off a series of independently chosen landing times, This series can then 
be taken as the actual landing times of successive aircraft for a hypothe 
tical aircraft carrier. Data for a large number of hypothetical landing 
patterns can be obtained in this way, and the empirical distributions, 
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means, standard deviations, etc., for various alternative landing sys- 
tems can be worked out. 

tAth this method of procedure it is clear that we have to decide 
on the basic frenuency distribution of landing times beforehand, though 
we could consider a limited range of types of distribution (rectangular, 
parabolic, etc.). 





In order to obtain a basic frecusncy distribution, observations must 
be made of actual landings with the landing times recorded. The usual 
curve fitting technicues can then be employed, Once the proper distri- 
bution has been decided upon (this may not be easy) the method of maxi- 
mum likelihood may be used to find the unimown parameters from the 
sample data (Mood 10 (۰ 

For computational ease, grouping of the frequency distribution could 
be utilized. From a practical point of view, a truncated distribution 
may be more realistic since the probability of extreme values are ef- 
fectly zero. After the preliminary planning, the corresponding dis- 
tribution function for each type of curve decided upon must be calcula- 
ted to four decimal places, Four figure numbers are then read from the 
table of random numbers, and the value of x (landing time) corresponding 
to the mid-points of the intervals into which the numbers fell are then 
recorded.. In this way series of landing times can be obtained readily, 
simply by taking successive groups of four figure numbers from the 
tables of random numbers and writing down the corresponding value of 
x read from the calculated table of the appropriate distribution function. 

Now for any series of landing times the waiting time of each air- 
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craft and the carrier's idle or wasted time at any time my be calcu- 
lated by the usual set of recursive equations (see Fig. 1.1): 


Tobi So V m .. 


If T,,, is negative, that is simply the carrier's idle tim, "tha 

constant arrival interval, the computations are simple but tedious. 
Finding the distribution of the actual stack is a little more com 

plicated. The easiest method of computation seems to be to consider 

for how long in the intervals (assumed constant) between successive 

arrivals there are 0, 1, 2,....,n aircraft waiting. For example, con- 

sider a constant departure interval of l unit., Consider five successive 

aircraft who have waits of 0, 0.5, 1.5, 2, and 1 units, Assume there 

is no stack when the first aircraft enters the landing pattern. Thus 

the first aircraft doesn't wait, so there are no aircraft in the stack 

and 1t is empty for l unit, say. The second aircraft waits 0.5 units 

of the second interval, thus the stack is empty for the remining 0.5 

units, The third aircraft waits sll of the third interval and 0,5 units 

of the fourth interval; while the fourth aircraft waits all of the fourth 

and fifth interval; and the fifth aircraft waits all of the fifth inter- 

val. Thus it follows that in the second interval there is one aircraft 

in the stack for 1 unit. In the fourth interval there are two then one 

alrcraft waiting for 0.5 units. It is not too difficult to carry out 

the arithmetical operations by using tables for visual representation 


as follows: 
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0.5 

iB 1 0.5 

2.0 0.5 1 

1.0 

The waiting times are grouped according to the arrival interval, The 
times during which O, l, ...., n aircraft are waiting during each inter- 
val are shown in the right hand table, and are obtained from successive 
diagonals of the center array. The figure appearing at the N, E, point 
of any diagonal shows how long the stack is occupied by the greatest 
mumber of aircraft in the corresponding interval, The latter number is 
given by the number of „uantities in the diagonal. The difference be- 
tween this figure and the one immediately to its S. W. gives the length 
of time that one fewer aircraft are waiting, etc. Subtracting the figure 


appearing at the S. y. point of any diagonal from 1.0 gives the length 


of time that the stack is wmoccupled. 
The above discussion explains how to obtain series of aircraft wait- 





ing times, carrier idle time, and the stack size. For any particular 
hypothetical carrier landing process one can calculate the mean aircraft 
waiting time; the total carrier idle time; and the tim distribution of 
the stack, By considering a series of such processes it is then possible 
to estimate what can be expected to happen in the long run, provided the 


observational landing data are valid. 
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